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We study the effects of instrumental systematics on the estimation of primordial non-Gaussianity using the 
cosmic microwave background (CMB) bispectrum from both the temperature and the polarization anisotropics. 
For temperature systematics we consider gain fluctuation and beam distortions. For polarization we consider 
effects related to known instrumental systematics: calibration, pixel rotation, differential gain, pointing, and 
ellipticity of the intrument beam. We consider these effects at the next to leading order, which we refer to as 
^— «j ' non-linear systematic effects. We find that if the instrumental response is linearly proportional to the received 

^v^j CMB intensity, then only the shape of the primordial CMB bispectrum, if there is any, will be distorted. We 

show that the nonlinear response of the instrument can in general result in spurious non-Gaussian features on 
• both the CMB temperature and polarization anisotropies, even if the primordial CMB is completely Gaussian. 

We determine the level for both the linear and non-linear systematics parameters for which they would cause no 
significant degradation of our ability to constrain the primordial non-Gaussianity amplitude /nl. We find that 
the non-linear systematics are potentially bigger worry for extracting the primordial non-Gaussianity than the 
linear systematics. Especially because the current and near future CMB probes are optimized for CMB power- 
spectrum measurements which are not particularly sensitive to the non-linear instrument response. We find that 
if instrumental non-linearities are not controlled by dedicated calibration, the effective local non-Gaussianity can 
be as large as /nl ~ O(10) before the corresponding non-linearities show up in the CMB dipole measurements. 
\ The higher order multipoles are even less sensitive to instrumental non-linearities. 
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C$ ; I. INTRODUCTION 

Characterizing the non-Gaussianity in the primordial perturbations has emerged as a powerful probe of the early universe. 
J> 1 Different inflationary models predict different non-Gaussian signal in both amplitude and shape (see review |Q]]). The specific 
C | departures from Gaussianity are highly model-dependent and any detection would rule out the simplest single field inflationary 
■ models. Non-Gaussianities measures the strength of interaction during inflation and can be studied and characterized most di- 
rectly by using the effective theory approach Jj-Hl- The first statistic that captures non-Gaussianity is the three-point correlation 
function or bispectrum in Fourier space. The overall amplitude of non-Gaussianity constrained from the data is often quoted in 
terms of a dimensionless non-linearity parameter /nl, and a shape which specify the configurations of the wavevectors of pertur- 
bations that contain the highest contributions to the non-Gaussian signal. With the assumptions of translational, and rotational 
invariance, the bispectrum shapes can be characterized as different triangle shapes. The three most studied shapes are: the 'local 
shape' where the non-Gaussian signal is maximum on squeezed configurations {k\ <K hi, £3), the 'equilateral shape' where the 
bispectrum peaks mostly on equilateral triangles (k\ ~ Ic2 ~ £3), and the 'orthogonal shape' which peaks for both equilateral and 
flat-triangle configurations (see e.g. JH|^] for a review). 

Different models of inflation predict different levels of /nl, ranging from (9(1) to /nl ~ (9(100). Values above /nl ~ (9(100) 
have been ruled out by the WMAP data already. Non-Gaussianity from "classical" inflation models that are based on a slowly 
rolling scalar field is very small [7, 8]; however, a very large class of more general models with, e.g., multiple scalar fields, 
features in inflaton potential, non-adiabatic fluctuations, non-canonical kinetic terms, deviations from Bunch-Davies vacuum, 
among others generates significantly higher levels of non-Gaussianity. Equilateral shapes are generically a signature of nonstan- 
dard kinetic terms in the inflaton Lagrangian, as for example in DBI |@| and ghost inflation jioll . 

Detections of any type of primordial non-Gaussianity would have profound implications on our understanding of the early Uni- 
verse and it is therefore crucial to estimate all possible contaminations of non-Gaussianity in CMB. Any physical/instrumental 
systematics to the non-Gaussian measurement must be well understood and controlled. We have to ascertain that any detected 
non-Gaussian signal has a physical origin and is not spurious. Potential contaminants include residual foreground contami- 
nation and unresolved point sources, non-Gaussianity induc ed by s econd-order anisotropies, such as gravitational lensing, the 
Sunyaev-Zel'dovich effect, and non-uniform recombination 
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In this paper, a detailed study of the impact of instrumental systematics on the primordial bispectrum measurement from 
CMB temperature and polarization fields is provided. We first show how such systematics change the CMB bispectrum signal 
for both temperature and polarization fields. We show that apart from distorting the primordial bispectrum, non-linearities in 
the instrument can also generate spurious bispectrum even in the absence of primordial bispectrum. This shape distortion could 
modify the effective normalization of a /nl estimator or even confuse the various primordial shapes of the bispectra. For a 
reliable inference of the level of primordial non-Gaussianity from upcoming CMB data, and thereby probe of the early Universe 
physics, precise and realistic estimation of instrumental systematics are required. 

The remainder of this paper is organized as follows. In §Q]]we introduce the basics of CMB bispectrum calculation. In jlTTI 
we present the parametrization of different instrumental systematics and explain the relations to the real experimental effects. 
In jjlVl we detail the modeling of linear and non-linear calibration systematics of CMB temperature anisotropies, and show that 
the systematics would generate new bispectrum and distort existing bispectrum. In |V]we calculate the systematics effects on 
bispectrum involving CMB £-mode polarization. Our numerical results and discussion are given in jjVIl where we also calculate 
the bias due to the non-linear calibration, the induced statistical error in /ml parameters due to non-Gaussian statistics of the 
distorted CMB fields. We compare the systematics requirement from bispectrum measurement with the requirement derived from 
the CMB B-mode polarization detection in jjVHI We also discuss our findings in the context of ongoing PLANCK and other 
upcoming CMB experiment. We summarize and draw our conclusions in Willi All our instrumental systematics calculations 
are carried out within the Jones matrix formalism where in Appendix A we describe the alternative, Stokes-parameters-based 
parametrization of beam distortions. In appendix B we describe CMB bispectra obtained in the full-sky formalism to supplement 
our calculations which have been carried out under the approximation of flat sky. 



II. THE PRIMORDIAL NON-GAUSSIANITY IN THE CMB BISPECTRUM 



To characterize the non-Gaussianity one has to consider the higher order moments beyond the two-point function, which 
contains all the information for Gaussian perturbations. The 3-point function, which identically vanishes for Gaussian perturba- 
tions, contains information about non-Gaussianity ITU pi \6L fO - Ua, for details]. Non-Gaussianities from the early universe can 
be described by the 3-point correlation function of Bardeen's curvature perturbations, which can be simplified assuming 
translational symmetry to 

<0(k 1 )0(k 2 )<I)(k3)> = (2w) 3 tf 3 (k, + k 2 + k 3 )/ NL ■ F(k u k 2 , h) (1) 



where F(k\,k 2 ,k-i) describes the shape of the bispectrum in Fourier space while the amplitude of non-Gaussianity is captured 
by the dimensionless non-linearity parameter /nl- The shape function F{k\,k 2 ,ks) correlates perturbations with three wave- 
vectors and form a triangle in Fourier space. Depending on the physical mechanism responsible for the bispectrum, the shape 
of the 3-point function, F(k\,k 2 , £3) can be broadly classified into three classes lfl7l,[l8ll . The local, "squeezed," non-Gaussianity 
where F{k\,k 2 ,kj) is large for configurations where k\ «: k 2 ~ £3. Next, the "equilateral," non-Gaussianity where F{k\,k 2 ,k-i) 
is large for configurations where k\ as k 2 ~ £3. Recently, a new bispectrum template shape, an orthogonal shape, has been 
introduced 1 19] which characterizes the size of the signal (f^" ho ) and peaks at both equilateral and flat-triangle configurations. 
It is easy to show that the shape dependence for the local model, fi oc .(&i, k 2 , £3), takes the following form 



Fioc.(ki,k 2 ,h) =2A% 



,3-Cv 
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where the normalized power spectrum Ao is defined in terms of the tilt n s and the curvature power spectrum as P<t,(k) = 
A< 6 £~ 3+( " ! ~ 1) with Po(k) defined in terms of the Gaussian component alone, <0 G (ki)<E> c (k2)> = ^\^-n)P<b(ki). 

Equilateral forms of non-Gaussianity arise from models with non-canonical kinetic terms such as the DBI action I20I1 . ghost 
condensation iflQIl . or any other single-field models in which the scalar field acquires a low speed of sound I2II l22tl . Although 
the shapes predicted by different models are in this case not identical, it has been noted lfl7l l23ll that they are all very well- 
approximated by the function 



F eq .(k u k 2 ,k 3 ) = 6A% 
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The definition for the equilateral model follows from the local one since is defined in such a way that for equilateral 
configurations, F eq ,(k,k,k) = Fi oc ,(fe, k, k) and one obtains the same value for given f^. = /Jjf '. The local and equilateral 
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forms are nearly orthogonal to each other, which implies that both can be measured nearly independently. The orthogonal form 
of non-Gaussianity is nearly orthogonal to both the local and equilateral forms lfl9ll 

F OIt ho(ki,k2,k?,) = 6Aj 



t 4 ti s 7,4-« T j4— n s iA—n s lA—n s iA~ n * 
K l K 2 K 2 K 3 K 3 K \ 



.(4-n s )/3,2(4-n.«)/3,4-« ! 

K 1 K 2 K 3 



+ (5 perm.) 



(4) 



The motivation of this shape is that a certain inflationary models yield a distinct shape which is orthogonal to both the local 
and equilateral forms. The orthogonal form has a positive peak at the equilateral configuration, and a negative valley along the 
elongated configurations 12411 . 

All single-field inflation models predict small amplitude for bispectrum in the limit of squeezed configurations, f^"' 1 = 
^(1 —n s ), regardless of the form of potential or kinetic term. Given that 1 -n s ss 0.04, all single field models predict f{%f al « 0.02. 
Hence the detection of local-type non-Gaussianity with amplitude greater than (9(0.1) would essentially rule out all single-field 
inflation models. For a credible non-Gaussianity detection, it would be essential to control all effects which might potentially 
contaminate the bispectrum. 

Primordial non-Gaussianity in the CMB: The harmonic coefficients of the CMB anisotropy af for temperature or Zs-mode 
polarization can be related to the primordial Bardeen curvature fluctuation <l>(k) via 



a x im = An{-i? § |^<D(k) «?(*)!£,(£), (5) 



where g x (r) is the radiation transfer function of temperature (X = T) or polarization (X = E). The most promising way of 
measuring non-Gaussianities is to study the CMB angular bispectrum, defined as follows 

r>XYZ _ ,X Y Z v 

D eie 2 (3,mim 2 m 3 ~ ^ a l\m\ u l 2 m 2 u e 1 ,m 1 ' ■ W 

Using Eq.©, the angular-averaged bispectrum can be written as 

a hhh ~ ( ^ n)( l > Zj m, m 2 m 3 J I (2tt) 3 (2tt) 3 (2tt) 3 ( l) { ™ K 2 ' f,m,( 3j 



x4(^i)4fe)4fe) «D(ki)<D(k2)*(k 3 )> (7) 

where (<t>(ki)<t>(k2)<D(k3)) is the primordial curvature three-point function as defined in Eq. (Q~|) and the matrix is the Wigner 3J 
symbol imposing selection rules which makes bispectrum zero unless 

(i) i i + T? 2 + ^3 =integer 

(ii) mi + ni2 + m 3 — 

(iii) \li - £j\ < 4 < ti + tj for i, j,k= 1,2,3. 

For simplicity, it is customary to define the reduced bispectrum b^^ z ( 



B XYZ _ y / i\ h h \gXYZ /o-v 



(2£i + l)(2€ 2 + l)(2tf 3 +1) ( £1 € 2 £3\,XYZ rQ x 



An \ 

Following convention, we define the radial functions a(r),B(r), 5(r), and y(r) for a given co-moving distance r for the temperature 
and Zs-mode polarization as 

of (r) = lfikk 2 g x e {k) j t {kr) , B x (r) = § ^fc *~ 1 g x (k) ;>(/tr)A , 

if (r) = § ^ MkrWf k(»>-W , d?(r) = \ fikg x (k) j e (kr)Af k^'\ {W) 

In the expression je(kr) is the Bessel function of order (. The transfer functions g T ( (k) and g^ (k) are numerically calculated using 
publicly available codes such as CMB fast JH and CAMB JH]. 
The bispectrum for the local case can be simply written as 

^[-^(^(^(^ + 2 perm.] . (11) 
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For equilateral template, 



= 6 J ^dr\-c$y)f t y)f t y) + 2 perm. +^(r)y[ 2 (r)^(r) + 5 perm. - 2^(r)^(r)^(r)] . 
Finally, the bispectrum for orthogonal shape can be cast in the form lfl9ll 

h fL = 18 f' ' 2dr -^S r )PlS r ^lM^ + 2 P erm - +^('")rJ 2 W^W + 5 perm. - |^(r)tf^(r)^(r) 



(12) 



(13) 



f7af Sky Limit: For simplicity we adopt the fiat sky approximation throughout. The anisotropies can be decomposed in spin-0 
and spin +2 plane waves in Fourier space 



ATXx) = f -^« r (ky k " 

[Q±iU](x) = J -^[a E (k)±ia B (k)]e 2i ^e^\ 



The Fourier transform of 3-point function is 

{a x (k l )a Y (k 2 )a z (k 3 )) = (27i) 2 B XYZ (k l ,k 2 ,k 3 )6 2 (k l + k 2 + k 3 ) 
and the bispectrum B{k\,k2, k 3 ) simplifies to 

B XYZ (k u k 2 ,k 3 )*b x ™ 3 , 

where |k,-| = ti. 

The bispectrum signal-to-noise is 

S_ _J_ 

N ~ 



(14) 



(15) 



(16) 



where the Fisher matrix T is given by Il27l - l30ll : 



v~i v~i C d i\ d t\ d £\ xyz i /^-\^ x> I /^-\^ Y> i /^-\^ z ' x'y'z' 

= xUikJ W v*) 1 W tAtA il > [ - } [ e J Wi 



(17) 



(18) 



CV C TE 



The sum over XYZ indicate permutation of TEE and TTE, and (C ( ) XY are the XY elements of the matrix ^ ^,f E j . We 

assume all the cosmological parameters, except /ml, are known. Indices XYZ and X'Y'Z' run over all the eight possible ordered 
combinations of temperature and polarization given by TTT, TTE, TET, ETT, TEE, ETE, EET and EEE. The overlap 
between any two bispectrum shape templates can be represented by a 2-d cosine as 12811 



cos(Bi,B 2 ) = 



B\ ■ B2 



{B\ ■ Bi) 1/2 (fi 2 ■ fi 2 ) 1/2 
where the dot product is defined as 

Bi B 2= f t l \ t\ f^ i^.4f3Mi,^,^/(C f ,Q,C f ,). 
J (2nY (2n) 1 (2n) 1 



(19) 



(20) 



The maximum magnitude of 2d-cosine, 1, correspond to a complete overlap between the two templates. Two completely 
orthogonal shapes will give rise to zero 2d-cosine. 



III. INSTRUMENTAL SYSTEMATICS 



The detection and measurement of CMB temperature anisotropy have already provided compelling evidence that the primor- 
dial perturbations have been generated during an inflationary period in the very early Universe. The next challenge is to constrain 
inflationary models by means of precise measurement of CMB temperature and polarization anisotropies. As we have noted, 
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the bispectrum of CMB is a powerful probe of the early universe. However, the instrumental systematic effects on the CMB 
bispectrum are relatively unexplored. 

For temperature systematics we can introduce the linear calibration (gain fluctuation) parameter of receivers, a(n). The non- 
linear response of the detectors to observed CMB fields parameter is parametrized by b(h). Effects associated with optical 
imperfections, i.e. beam shape, are collectively captured by the parameter c 

& ohs (h) = [1 + a(n)]©(n) + £(n)©(n) 2 + [c ® ©](n) + ... , (21) 

where © ois (n) = AT(n)/T is observed CMB temperature fluctuation at direction n including detector noise. We will see the 
effect of these systematics on the temperature bispectrum and eventually on /nl in the next section. The last term accounts for 
effects associated with beam systematics (or more generally, any effect that can be described as a convolution with the underlying 
sky in real space). Here ® denotes 2D real-space convolution. One immediate consequence of this is that by the convolution 
theorem this will correspond to a product in multipole space, the last term then reads c(l)0(l) which implies no / - / coupling 
in multipole space. Note that this argument neglects the effect of scanning strategy, which we assume is small throughout. This 
is well-motivated for temperature anisotropy measurements. However, this is not necessarily the case when the (very weak) B- 
mode polarization is considered and beam imperfections (such as differential pointing or differential beamwidth) leak the much 
larger temperature anisotropy to B-mode polarization via a scanning-strategy-dependent local weight which enters the B-mode 
map-making process (as well as maps of T and Zs-mode) and generally varies across the sky and depends solely on the non- 
uniform scanning strategy. Since only little information is contained in the B-mode polarization on primordial non-Gaussianity 
we will ignore this effect here. In principle, leakage of temperature anisotropy to Zs-mode polarization cannot be ignored and 
indeed we give a special attention to these terms here. However, we ignore the impact of the coupling of non-uniform scanning 
strategy to beam distortions on the spurious Zs-mode polarization and its propagation to /nl bias as this is a higher order effect 
and highly depends on the details of scanning strategy. 

Polarization systematics: The measurement of CMB polarization is complicated by its low level signal, galactic foregrounds, 
and various systematic effects. In particular, precise control of systematics is required, especially for polarization measurements. 
Measuring the polarization using Polarization Sensitive Bolometer (PSB) pairs inside a feedhorn, each measuring the intensity 
of one direction of the polarization, is used by several ground- and balloon-based CMB experiments such as Planck 1 , MAXIMA 2 , 
Boomerang 3 , BICEP 4 , EBEX 5 etc. The difference of the received intensity from the two PSB gives a combination of Q or U 
in the frame of the focal plane. Polarization measurements therefore involve taking the difference of received intensity at two 
orthogonal polarizers. Any differences between the two PSBs might generate spurious Q and U signals. Furthermore, the beams 
of the horns generally have some degree of ellipticity. The typical differential ellipticity of a pair is around few percent, while 
for a single horn is usually at sub-percent level ~ 0.5% ifJlll . When several detectors are combined to obtain Q and U maps, it 
is important to match the responses of these detectors precisely in terms of cross-calibration, beam shape, spectral response, etc. 
Below we discuss the parametrization of CMB polarization systematics. 

The polarization systematics fall into two categories, one associated with the detector system which distorts the polarization 
state of the incoming polarized signal (Type I hereafter), and another associated with systematics of the CMB signal due to 
the beam anisotropy (Type II hereafter). Each type of systematics contamination can be further classified into linear and non- 
linear effects. We parametrize the linear instrumental systematics for CMB polarization measurements following ll32ll (HHZ 
from hereafter), and introduce corresponding analogous non-linear systematics parametrization. In Appendix A we briefly 
discuss an alternative parametrization employing a Stokes-parameter-based formalism and show the equivalence of the two 
parametrizations. 

The instrumental response to incoming CMB radiation is usually described by the Jones transfer matrix. Bias induced in the 
matrix determination will mix the corresponding Stokes parameters. To first order, the effect of Type I systematics on the Stokes 
parameters can be written as ll32tl 

6[Q ± iU](h) = [a± i2co](n)[Q ± >U](h) + [fi ± ifiW)[Q + iU](h) + [ n ± iy 2 ](n)&(A) , (22) 

where a(n) is a scalar field which describes the miss-calibration of the polarization measurements, w(n) is another scalar field 
that describes the rotation miss-alignment of the instrument, (f\ + //iXn) are spin ±4 fields that describe the coupling between 
two spin +2 states (spin-flip), and (y\ + /y2)(A) are spin +2 fields that describe monopole leakage, i.e. leakage from temperature 
anisotropy to polarization. 



1 http : //www . rssd . esa . int/index . php?proj ect=Planck 

2 http : //cosmology . berkeley . edu/group/cmb/ 

3 http : //cmb . phys . cwru . edu/boomerang/ 

4 http://bicep.caltech.edu/ 

5 http : //groups . physics . umn . edu/cosmology/ebex/ 
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Similar to the Type I systematics, the effect of Type II systematics on the Stokes parameters can be written as [3J 
6[Q + iU](h; cr) = <xp(n) • V[g + z'E/](fi; cr) + a[d x ± id 2 ](S)[di ± ;d 2 ]0(n; cr) + <x 2 ^(n)[3i ± id 2 ] 2 0(n; cr) , 



(23) 



where the systematic fields are smoothed over the average beam cr of the experiment. Therefore, the type II systematic fields are 
sensitive to the imperfection of the beam on the beam scale, cr. The spin + 1 fields, (p\+ip 2 ) and {d\+id 2 ), describe pointing errors 
and dipole leakage from temperature to polarization, respectively, and q is a scalar field that represents quadrupole leakage l32il . 
e.g. beam ellipticity. 

We parametrize the non-linearities of the instrument as 



5[Q + iU](n) = [y x + if 2 ](h)®(h) 2 + cr[di + id 2 ](n)[d 1 + /<9 2 ]0 2 (n; cr) + o~q(h)[d x ± /<9 2 ] 2 2 (n; cr) 



(24) 



where j\ and y 2 are spin +2 non-linear parameters which describe the second-order leakage from temperature to polarization 
states; d\ and d 2 are spin +1 parameters which describe the dipole leakage from second order temperature to polarization states. 
q is a scalar field that represents quadrupole leakage from second order temperature. In general, we use S for systematics of 
non-linear response to distinguish the same type S linear systematics. Note that in principle there can be terms proportional to 
(Q + iU) 2 (h), etc, which we have dropped as the 2 (n) terms provide the dominant non-linear contribution. 



IV. SYSTEMATICS-INDUCED CMB TEMPERATURE BISPECTRUM 



In this section, we discuss the impact of instrumental systematics on the measured bispectrum. We show that the linear 
calibration systematics of detectors can only distort the primordial CMB bispectrum; any detection with only the presence 
of linear systematics implies the detection of primordial non-Gaussianity. However, non-linear systematics can generate new 
spurious bispectrum even in the absence of vanishing underlying bispectrum, i.e. the primordial CMB is purely Gaussian. As 
discussed before, beam shape effects will only modify the /'-dependence of an already existing bispectrum, i.e. in case that 
/nl = these beam effects are irrelevant. 

With both linear a(h) and non-linear b(h) systematics included, as in Eq. d2Tl ). the CMB 3-point function of the temperature 
can be written as 



<0 oii (n 1 )0 o6s (n 2 )0 o " s (n 3 )) CMfi , jy 



<0(ni)0(n 2 )0(n 3 ))cMB + <a(ni)> JyJf <0(ni)0(n 2 )0(n 3 )>cMB + per m 
+<Kni)>™,<0 2 (ni)0(n 2 )0(n 3 )>cMB + perm 
+<fl(ni)a(n 2 )> JVJ ,<0(ni)0(n 2 )0(n 3 ))cMB + perm 
+<a(ni)fe(n 2 )). !J ,. ! ,<0(ni)0 2 (n 2 )0(n 3 ))cMB + perm 
+<a(ni)a(n 2 )a(n 3 )) jyj ,<0(ni)0(n 2 )0(n 3 )>cMB 
+</7(ni)fe(n 2 )). SVJ ,<0 2 (ni)0 2 (n 2 )0(n 3 )>CMB + perm, 



(25) 



where {...)cmb means averaging over CMB realizations, and (...) svsf stands for averaging over systematics field realizations. 
In principle, there can be higher order terms proportional to 3 (n). For 3-point correlation functions, we are only interested 
in terms up to second order in CMB fields. However, for the effects on trispectrum estimations, 3rd order terms should be 
considered. The first term on the right hand side of Eq. d25l l is the cosmological 3-point correlation function, i.e. the target 
signal. The second term is proportional to linear systematics field, e.g. first order gain response. If there is non-zero monopole 
gain systematics, (a(n)) syst , then this term will induce a linear bias proportional to the primordial bispectrum. For the third 
term, if there is any non-zero (b(ik)) syst , then this term will introduce a new bispectrum even if there is no primordial non- 
Gaussianity. The forth term contains self correlations of linear distortion field parameter and is capable of only distorting 
existing non-Gaussianity. In contrast, the fifth term coming from cross-correlation between linear and non-linear systematics, 
can generates new bispectrum. If the systematics field is non-Gaussian then the term containing three distortion parameters (the 
term proportional to (a(ni)a(n 2 )a(n 3 )) svsf ) can also contaminate /nl measurements. We do not consider this term and terms 
involving (b(hi)b(n 2 )) „,„ in this paper. 

In case that a sufficiently small sky patch is considered, spherical harmonic modes can be replaced by Fourier modes. Gener- 
alization from the flat-sky to the full sky is straightforward and is given in Appendix B. If we consider the non-linear calibration 
up to second order, the Fourier transform of the Taylor-expanded observed temperature anisotropy (including calibration con- 
tamination) could be written as 



&(l) ohs = f rfn0*(nK M = 0(l)- f -^r0(l') 

J J (2n) z 



r d 2 i" 

-»')+ ^0d>(l-l'-l") + 

J (2tt) z 



Using the flat-sky approximation, we can define the angular bispectrum as 

<0(1 1 )0(1 2 )0(1 3 )) = (2^) 2 C T(1 1 + 1 2 + 1 3 )< 5 3) . 



(26) 



(27) 
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Linear Systematics S 


W s E (\ u h) 


Calibration a 


cos[2(^i, - ^l)] 


Rotation a> 


-2 sin[2(<,oi 2 - ipjj] 


Pointing p a 


cr(h ■\ l )sm[2((p h _ - <p L )] 


Pointing pt 


-cr(l 2 X Ii) • t sin[2(vi 2 - <Pl,)] 


Flip fa 


cos[2(2v>i, - (fi\ 2 - <Pl)] 


Flip 


- sm[2(2tp h - <pi 2 - tp L )] 


Monopole y a 


cos[2(<ft, - tpi,)] 


Monopole y b 


- sin[2(<p h - L ) 


Dipole d a 


(l z &) sinfrffi, + </> h - 2ip L ] 


Dipole dt 


(l 2 cr)c(K[tp h +(p u _ - 2<p L ] 


Quadrupole q 


-(hcrf cos[2(<£i 2 - </> L )] 



TABLE I: Window functions for all the 1 1 linear response systematic parameters. First column indicates the type of systematic parameter in 
consideration. Second and third columns show the window functions needed to calculate systematic effects on bispectrum measurement and 
appear in Eq. d33t . We note that li = /Ji, h = L — li, and li = hh and ip\ = cos _l (fi • 1). 



The effect of systematics on bispectrum can be obtained by Fourier transforming Eq. d25l ). For simplicity we will classify 
the effects from systematics into two categories; one where the systematics create new bispectrum even if the underlying CMB 
is purely Gaussian, referred to as a 'bias term', and second, where the shape of the underlying bispectrum is distorted due to 
systematics. The first contribution can be written as 



B 



QQQMas 
(JiJbJb) 



r a@ r @@ 

W2 



5 perm. 



b(0)Cf & Cf & + 5 perm. 



r d 2 v 

J (2nf 



C] _yC f , C t + 5 perm. 



-,00^00 



(28) 



The first term is zero since one does not expect systematics to be correlated with the CMB. This term is analogous to the 
integrated Sachs-Wolfe (ISW) effect which cross correlate the lensing potential with CMB field. The second term is proportional 
to b(0) which is the bispectrum generated from the monopole of non-linear systematics field. The terms involving the cross 
correlation between the a(l) and b{\) fields can induce a bispectrum as well. We note that with no cross-correlation between the 
linear a{\) and non-linear b{\) field, only the monopole of the non-linear parameter b(0) would bias the bispectrum detection. 
The term describing the systematic distortion of primordial bispectrum can be written as 
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The total observed bispectrum B?® '^ would then be 
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(29) 



(30) 



V. SYSTEMATICS-INDUCED CMB POLARIZATION BISPECTRUM 

Currently the best CMB constraints on the three parameters f^[ al , ff"' 1 ', come from the WMAP temperature anisotropy 
data [fl9L L33l - l36h . By also having the polarization information, one can improve sensitivity to primordial fluctuations ITU Eii 
I29I l37ll . The ongoing Planck experiment and futuristic CMB experiments such as CMBPol will characterize the polarization 
anisotropy to high accuracy. In this section we discuss the instrumental polarization systematics which can contaminate the 
primordial non-Gaussian signal. 



Impact of linear response on polarization bispectrum 

Here we show the calculation for spurious bispectrum distorted by linear instrumental systematics involving polarization 
field. Performing the harmonic transformation of Eq. (f22b and Eq. (l23l l. we get the distorted polarization E'-mode from linear 
instrumental systematics 

E{\) obs = j dn [Q obs (a) cos(2^i) - U obs (h) sin(2^)]e"' l fi . 

(31) 
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Non-linear Systematics § 


W£(l,l',l") 


Monopole leakage y a 
Monopole leakage % 
Dipole leakage d a 
Dipole leakage dt 
Quadrupole leakage q 


2 cos[2(^m-_i» - <pi)] 
-2sin[2(^,_ F _ r - ip{)] 
cr[l' sin(^, + ifii-v-v ~ 2<P\) + 1" sin(^y + <py- V -\» - 2^)] 
cr[l' cosftcy + ^_ F _!« - 2<p x ) + 1" cos(^,« + - 2<^])] 
-2o- 2 jl' 2 cos[2( Vl , - + l" 2 cos[2( Vl „ - </>,)]) 



TABLE II: Window functions for 5 non-linear response systematic parameters. Here li = /Ji, L = L - li, and I2 = hh, and ip\ = cos '(n • 1). 



We can easily calculate the bispectrum involving observed Zs-mode to find out the changes induced by different types of linear 
instrumental systematics 

(X obs (h )Y° h \h)Z ob \h)) se (2n) 2 6(h +h + h)B^l w (32) 

where X, Y, Z is temperature or E-mode polarization field of CMB. For example, the bispectrum SB^f^f describing the distor- 
tion of primordial polarization bispectrum due to instrumental systematics can be written as 

+B^v h+ vM)W S E (}n -l')W|;a 3 , -1") + B^ lM _ lr) W s E &, -Y)W s E ,(h, -1")] . (33) 

Summations are over 1 1 systematics contributions, S = a,a>, f a , fb, p„, pb,y a ,7b,d a ,db,q- Window functions I2) are given 

TableU 



Impact of non-linear response on polarization bispectrum 

As for the linear response systematics, we can calculate the observed CMB £ -polarization in the presence of non-linear 
systematic fields. We focus on the non-linear leakage from second order temperature to £ -polarization, as this is the dominant 
systematics which could potentially bias bispectrum detection. We will show that the monopole of the non-linear response 
systematics 71,72,^1,^2, and q fields can generate new bispectrum. To second order, any cross-correlation between linear and 
non-linear systematics could also generate non-zero bispectrum. The £(1) field in the presence of the non-linear systematics § 
can be simplified as follows 

E° b \Y) = E obs (\) + f f ' f 5(1 - 1' - l")f (fit (1")W| (U',1") • 

Unlike the temperature bispectrum case, for polarization there are many terms containing cross-correlation between non-linear 
systematics (7, d, or q) and linear systematics S (11 systematics as given in Table [Q. The spurious bispectrum contribution for 
EEE, TTE, and the TEE bispectrum can be written as 

B ZhM laS = S (W C l Ec l E W l 03 . "li. -I2) + Perm. 
Klhif aS = S (OX^C™" w£ C»3. -»i • -!2> + P erm - 

B Zkl) iaS = 5(0)C f 7 ; r C[ 2 £ w|(l3,-li,-l2) + 5(0)C, 7 ; r C[, £ w|(l2,-li,-l3) + perm. (34) 

where W E are given in Table [TT] The total observed bispectrum B^^°J including both the systematics induced and distorted 
contribution is 



r>Alz.,oDs D XYZ,primordial ST 1 XD XYZ,dist . ST 1 nXYZMas 

B - B (l„l 2 ,l3) + h 5B (h,h,h) + Zj B (h,hM ■ (35) 



nXYZ,obs 
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VI. NUMERICAL RESULTS 



In order to numerically calculate the systematic effects on bispectrum, one needs to choose specific model of systematics 
field, i.e. the power spectrum C s e s of the systematics field S . So far we have not made any assumption about the systematic 
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FIG. 1 : Upper left panel: Spurious bispectrum generated from the cross correlation of calibration parameter a(L) and non-linear response 
parameter b(L). For reference the solid red curve shows the primordial bispectrum of local type with non-linearity parameter /nl = 1. Upper 
right panel: Distortions on temperature local type bispectrum from linear instrumental gain fluctuations. For reference the solid red line 
shows the CMB bispectrum from local type model with / NL = 1. Lower left panel: Similar for upper right panel but for equilateral type 
of bispectrum. Lower right panel: Similar for upper right panel but for orthogonal type of bispectrum. In all the panels We have assumed 
systematics model as defined in Eq. i37\ and choose the rms of systematics field A a = A b = 0.5, while varying coherence length starting from 
o~s = 10' to o~s = 120'. Note that the spurious signals scale as A 2 S . 



field. In reality the power spectra of the systematics fields will depend on the experimental design and scan strategy, and could 
be inhomogeneous, anisotropic and complicated. However, we here employ a model which, although not exact, can be used to 
assess the level of contamination, independent of the origin of these systematics fields. As a simple model we assume that the 
contamination fields are statistically isotropic and Gaussian, thus their statistical properties can be fully described by their power 
spectra ll32tl . 

(S (1)S '(1')) = (2tt) 2 (5(1 + l')C s e y , (36) 
where S stands for any of the systematic fields. We assume the power spectra of the form 

c ss = A l expH0 + l)gf) 

i.e. white noise above certain coherence scale cr s , which is a key quantity that sets the level of contamination of each systematics 
effect together with As , which describes the rms of the contamination field S . Note that <x without the subscript S refers to the 
instrumental beam FWHM. 

In Fig ([]]) we show the systematics-induced temperature bispectrum B^^f whose analytic result are given by Eq. ( l28b 
and j29l . In the upper right and lower panels we show the effect of linear systematics a which distort the CMB bispectrum of 
local, equilateral, and orthogonal templates. In each case we consider several choices of coherence length varying from cr s = 10' 
to cr s = 120'. In the upper left panel we also show the systematic-induced bispectrum from the cross term C" b of Eq. d28l ). 
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FIG. 2: Upper panels: The requirement on non-linear systematics parameter b at lcr level as a function of FWHM for various choices of 
instrumental sensitivity and for local (left), equilateral (middle), and orthogonal (right) templates. The requirement were derived by demanding 
that the spurious local f£T' generated by systematics are detectable at lcr i.e. by solving Eq. ( I39l >, Lower panels: The minimum detectable 
effective f^' from non-linear systematics for b = 1 as a function of FWHM and for 3 choices of instrumental sensitivity, for local (left), 
equilateral (middle), and orthogonal (right) templates. 



The systematics-induced bispectrum scales as A,. In all the these panels the nns amplitude (A a — At — 0.5) was set to high 
values for illustration purposes. For experiments with systematics nns fluctuations As ~ 10%, the linear distortion is not the 
dominant effect for i < 2000 if coherence cr s > 10'. For experiments which probe much smaller scales (£ > 2000), the extracted 
non-Gaussianity might be contaminated by linear systematics. Note that for large coherence lengths both distortion bispectrum 
and bispectrum induced due to Of have a similar shape as the primordial bispectrum. However, for small coherence length the 
convolution transfers power from large to small scales and the resulting distorted bispectrum is much flatter than the primordial 
bispectrum. For large coherence length the systematics power spectrum drops rapidly and does not have power over broad 
enough I range to transfer power around. 

We want to quantify the maximum tolerable systematics (for a given experiment) below which they would cause no significant 
degradation of our ability to constrain primordial non-Gaussianity. We again focus on three templates, local, equilateral and 
orthogonal. The confusion of non-Gaussianity /nl from systematics with the primordial non-Gaussianity can be quantified by 
solving 

Z V f ffil ^ 2 ^ 3 DXYZ,syst(^-l\ XX ' (^-l\ YY ' (^-l\ ZZ 'x'Y'Z',prim 

f syst V V f d t\ d h d lj XYZ ,prmi(^-l\ XX ' (^-l\ YY ' (^-l\ ZZ 'x'Y'Z\p f wh-, 

/nl^J(^^(^ b w3 (QJ (QJ (QJ K** W> 

where C xx ' is the observed CMB power spectrum with noise. Here f®£ is an effective amplitude of non-Gaussianity generated 
for a given template shape. We calculate the requirement on non-linear systematics errors by demanding that the effective fj^' 
is not detectable at lcr using the optimal CMB estimator as developed in Ref. I23ll29l[35l[38ll . i.e. jS£' < (/nl 2 ) 1 ^ 2 , giving the 
requirement 



I: 



(nshape,prim\ 
B V2'3 ) 



- qWJP^qW Ci.Ct.Ce, 
b shape < _J ' 2 1 . (39) 

n/uVi r nnniu.s ■ nave x ' 



d 2 h d 2 l 2 d 2 h 

(2m) 2 (2k) 2 (In) 1 Ct.Ct.Ct, 



In Fig. (f2]i upper panels we show the maximum allowed non-linear temperature systematics b as a function of beam FWHM for 
local, equilateral and orthogonal shapes and for various choices of experimental noise sensitivity. For a cosmic-variance-limited 
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FIG. 3: 2-d cosine between the non-linear systematics induced bispectrum and local, equilateral, and orthogonal template. The left panel is 
for the temperature systematics bispectrum B TTT - b and the middle and right panels are for the polarization B EEE -' 1 , and B TTE &, respectively. For 
the polarization plots (middle and right), the spikes in the plots correspond to sign-flip in the cross-correlation of systematics and primordial 
bispectrum. 

experiment up-to i max = 2000, the requirement for local, equilateral and orthogonal shapes are b < 22 (2cr), b < 41 (2<x) and 
b < 44 (2cr) respectively. We will see in the next section that these requirements from bispectrum are more stringent than the 
requirements for non-linear systematics from the fi-mode power spectrum. 

In Fig. (O lower panels show the effective fS£' from the non-linear systematics for local, equilateral and orthogonal templates 
as a function of beam FWHM and for several choices of experimental noise sensitivity. Note that systematics-induced f^' is 
only weakly sensitive to the instrument beam and detector sensitivity. Also since f5?' is of order unity for b — 1 , the requirement 
on b is comparable to minimum detectable /nl i.e. /nl ~ b for given CMB experiments. As expected the requirement on 
the systematic control gets weaker as the FWHM and/or instrumental noise is increased. Also note that since the maximum 
multipole used in the analysis was fixed to C max = 2000, all experiments with low enough instrument noise and FWHM give 
identical requirement because they all are cosmic variance limited up-to C max = 2000, as indicated in Fig. (0. 

To quantify the correlation between the systematics-induced bispectrum and primordial bispectrum template, one can calcu- 
late the 2-d cosine as defined in Eq. ( fT9b . Fig. ([3]) shows this 2-d cosine, between the systematics and local, equilateral and 
orthogonal templates. Systematics have comparable overlap with all the three templates, and overlap is relatively small with 2-d 
cosine of order ~ (9(0.2) for equilateral and orthogonal shape at i max ~ 2000, and of order ~ (9(0.4) for local shape. Since we 
modeled the temperature non-linear response as 0(n) = @(n) + M9(n) 2 , one might expect a large overlap with local template 
however this not the case because (1) the local template is defined for three dimensional k space of primordial curvature pertur- 
bations, while we show the overlap with the projected two dimensional ^-space bispectrum; (2) at high t there are oscillatory 
features in both the template and systematics bispectrum which after summation over I cancels out. Because of the week over- 
lap between systematics and primordial templates, systematics generate comparable amount of non-Gaussianity for the three 
templates considered. 

Given the systematics contamination ' for a given shape, one can get a rough estimate for the level of f5?' for another 
shape by using the 2d-cosine between the two shapes. For example, if we start from systematics contamination for equilateral 
shape for a cosmic variance limited experiment, f^' equl ~ 2.6 (see Fig. (2) lower middle panel). We can convert it to the 

effective local f^'' hcal by using a fudge factor ~ 6 between the two shapes to get / N "'' /ocfl/ ~ 2.6/6 = 0.44. This agrees nicely 
(see Fig. 2 lower left panel) with what we get with the full calculation. The fudge factor 6 comes from 2D cosine which gives 
a factor of 2 and the normalization which gives a factor 3 fl7ll . The reason we get higher /nl in the equilateral model is that 
we normalize the maximum of this shape to the minimum of the local shape. The normalization factor 3 compensates for this 
effect 6 (see Ref. JH] for details). 

For polarization, many of the instrumental systematics can be removed or suppressed by carefully designing the scan strategy. 
Instrumental rotation is potentially a powerful way to mitigate instrumental systematics effects if the systematics field has 
different spin-dependence than the CMB polarization Q ± ill. Ideally, if every pixel can be covered by infinite number of 
different rotation angles with respect to the real polarization direction on the sky, systematics fields (y\,y2,d\,di,u),f\,f2,p\, 
and P2) which have different spin from Q + iU can be completely removed. However the quadrupole leakage (leakage from 
temperature to polarization, stemming from beam ellipticity) has the same spin dependence as the polarization parameters 
Q + iU and hence can not be averaged over, even with the perfect coverage of rotation angles. For this reason we mainly focus of 
quadrupole leakage q. As with the temperature systematics b, we find that only the monopole of systematics field can generate 
new bispectrum. Fluctuations can only distort the existing bispectrum without generating spurious signals. 



The local and equilateral shapes are normalized to the same value at equilateral configuration, i.e. F/ oc „;(k, k, k) = F equ n(V., k, k) 
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FIG. 4: Effective / NL of local (left), equilateral (middle), and orthogonal (right) templates generated from non-linear quadrupole leakage q as 
function of FWHM (upper panels) and l max {lower panels) for EEE, TEE,TTE estimators. For lower panels we assume a cosmic-variance- 
limited experiment. The spikes in the plots correspond to sign-flip in the cross-correlation of systematics and primordial bispectrum. 



In Fig © we show the effective non-linearity parameter f^' generated by quadrupole leakage q for the local, equilateral and 
orthogonal templates. We show the fj?* ' as function of FWHM (upper panels) and as a function of maximum multipole used 
in the analysis t max . Since the quadrupole leakage come from the beam effects, larger the beam FWHM, larger the spurious 
non-Gaussianity fS['. The spikes in the /nl vs t max plots correspond to sign-flip in the cross-correlation of systematics and 
primordial bispectrum. 

In ([5]) we show requirements on non-linear polarization systematics parameter q using separately the TTE, TEE and TEE 
bispectrum. As for the temperature case, the requirements were derived by demanding that the spurious local f^[' generated by 
systematics are detectable at lcr. Note that from the definition of systematics leakage the requirements for quadrupole leakage 
q are expressed in the units of square of beam FWHM. In comparison to the temperature systematics the requirement from 
polarization is much weaker. Since the polarization systematics bispectrum contains sines and cosines, the dot product between 
the template bispectrum and polarization bispectrum oscillates, resulting in weaker effective f2['- Since the quadrupole leakage 
scales as q ~ a 2 , the effective shown in Fig dU also has the same scaling. 



VII. IMPLICATIONS FOR CMB EXPERIMENTS 



How well can non-linear systematics be controlled? The systematics requirements on the linear response parameter a for 
CMB experiments aimed at B-mode observations is well discussed in the literature ll32ll39ll 



< 0.06 { - j for coherence cr a — 10 

< 0.05 (—^— V /2 for coherence cr a = 120' (40) 

v 0.005 ; 



To estimate the level of non-linear response b which can be detected at the power-spectrum level, we start with 

6® obs (n) = (1 + a)6& sky (h) + b(6® sky (A)f , (41) 

and note that the non-linear term is negligible for the limits when 6& sky — > and contributes the most when S® sky = 5® max . 
Therefore, given the maximum allowed linear response parameter a, one can obtain a requirement on b as 
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FIG. 5: Requirements from the TTE, TEE and EEE bispectrum for the non-linear quadrupole leakage q as function of FWHM (upper panels) 
and C max (lower panels) for noise sensitivity A p = l^K-arcmin and for local (left), equilateral (middle) and orthogonal templates (right). 
The requirements were derived by demanding that the spurious local / NL generated by systematics are detectable at lcr. Note that from the 
definition of quadrupole leakage, the requirements are expressed in the units of a 2 and therefore the requirement becomes more stringent as 
the beam FWHM is increased. For lower panels the beam FWHM is 10', and the spikes correspond to sign-flip in the cross-correlation of 
systematics and primordial bispectrum. 



where A is a fudge factor of order unity. Primary CMB anisotropy will provide an upper bound of b « a x 10 5 from Eq. (|40T > l32ll . 
The requirements on a are more stringent from B-modes than from bispectrum systematics. Using the numbers we get 

I r W 2 

b < 6000 for coherence cr a — 10 

1 0.005 ; 

b < 5000 {q^) 12 for coherence cr a = 120'. (43) 

However, we can go beyond this rough estimate; one can use the dipole of the CMB to set more stringent requirement on b. 
With ©dipole ~ 0.001 one obtains a requirement which is ~ 100 times more stringent than provided by Eq. d43l . For r = 0.005, 
this requires b < 50, which is still less stringent than the requirement from bispectrum. Therefore, if instrumental non-linearities 
are not controlled, the effective non-Gaussianity can be as large as /nl ~ 30 (see figure 2) before the corresponding non- 
linearities show up in CMB dipole. The higher multipoles of the CMB power-spectrum are even less sensitive to experimental 
non-linearities. 



Non-linear response of the PLANCK satellite 

PLANCK is the third generation space radiometer, following the COBE-DMR and WMAP instruments, and is currently 
collecting data from the Lagrangian point L2. PLANCK is expected to significantly improve our knowledge on non-Gaussianity 
of the CMB. Here we discuss the the non-linear response of the Planck detectors, following Ref. ll40M42n . Two state-of-the- 
art and complementary instruments are integrated in the focal plane of the PLANCK 1.5-meter dual reflector telescope. The 
Low Frequency Instrument (LFI), based on coherent receivers, observes the sky in three bands centered at 30, 44 and 70 GHz. 
The HFI, using bolometers cooled down to 0.1 K, covers six channels between 100 and 850 GHz. The LFI is an array of 22 
cryogenic coherent differential radiometers based on indium phosphide HEMT (high electron mobility transistor) low noise 
amplifiers. The receiver array is split into a front-end unit (cooled down to 20 K to optimize sensitivity) and a back-end unit 
(operating at ~300 K). The instrumental properties of the LFI and HFI have been calibrated and tested at different integration 
levels by detectors, individual receivers, and the whole receiver array during ground calibration and in flight calibration. For 
LFI, the response of 30 and 40 GHz radiometers show slight output compression, which affects the noise properties. The effect 
on noise temperature, white noise sensitivity, and noise effective bandwidth has been studied extensively. Here we focus on the 
implications on bispectrum measurement. 
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The non-linearity of microwave receiver can be introduced by various sources, e.g. radio-frequency (RF) amplifiers, detector 
diode, and back-end analog electronics ll40ll . The back-end RF amplifiers and detector diodes might cause the non-linear response 
of 30 and 40 GHz receivers on LFI. For the LFI pseudo-correlation receivers, assuming the sky signal and the reference load can 
be perfectly isolated after the hybration 

^out G(T[ n , T^noise ) x (T iD + 

-^noise ), (44) 



where T m refers to any input signal, could be either T s k y or T le f, V out is the corresponding voltage output, r n oise is the corre- 
sponding noise temperature and G{T^ a , r no j se ) is the photometric calibration constant (detector gain) which may depend on the 
input and noise temperatures (non-linear response) in general. The parametrization of the non-linear response model has been 
provided in iPioll 

Vout = Gtot (7"in + TWse) > 

Gtot = 777777; — 77 : > (45) 

1 + kG (i in + r noise ) 

k is the photometric calibration constant, and have been estimated with various sky-load temperature for the 30 and 44 GHz 
receivers. If k = the response parameter (receiver gain), G tot , reduces to the constant Go, and k = 00 corresponds to infinite 
compression of the input signal. In this model, Go is proportional to gains of front-end module, band-end module, and the 
bandwidth. 

The linearity of the receiver response is important for the in-flight calibration. For example, photometric calibration is per- 
formed by exploiting the dipole anisotropy and for the beam measurements through observations of bright point sources, e.g. 
Jupiter and Saturn. 

The variation on G, SG/G, caused by a variation of the input temperature ST can be calculated as 



^ = k J£2l . (46) 

G l+kG (T sky + 6T + 

For dipole anisotropy, ST ~ +3 mK, and for bright point source like Jupiter, ST ~ ±50 mK. The estimation of the receiver 
parameters Go, ^noise and k can be found in l40ll . The non-linearity or the relative change on the detector gain has been estimated 
to be 



f <6x 10~ 5 for ST ~ +3 mK 



^ < 10~ 3 for ST ~ +50 mK, 



(47) 



With the above limits on the gain, we can now calculate the non-linear response of temperature field b of our parametrization 
by noting that SG/G = bT/(l + a + bT). For T ~ ±3 mK this gives b < 0.1. We can see that during nominal operation 
with small input signal, the performance of the receiver is sufficiently linear as not to contaminate the bispectrum measurement. 
For bolometer detectors like HFI, the dynamic range is much smaller. This is usually not a problem if low temperatures are 
considered, but for the simultaneous measurement of galaxy and the CMB the non-linearity parameter b > 0(10) might pose a 
problem. To our knowledge, there is no available literature which characterizes the non-linearities of the HFI instrument Il43m7ll . 
Since HFI will be the driving instrument for primordial bispectrum measurements, it is essential that instrumental non-linearities 
satisfy our benchmark criteria. 



VIII. SUMMARY 



Future CMB experiments will reach the raw sensitivity to detect local shape /nl ~ 5(l<x) using CMB temperature anisotropy 
information, and /nl ~ 2(1 cr) using combined the temperature and E'-polarization information. However to reach this Fisher 
limit and reliably infer the level of primordial non-Gaussianity, precise and realistic calculations of the effect instrumental 
systematics on bispectrum are required. Although the impact of instrumental systematics have been extensively studied for 
the CMB power-spectrum, the effect on bispectrum have not been extensively explored. We have studied, in a fairly general 
manner, the scientific impact of the instrumental systematic effects on the primordial non-Gaussianity measurement from the 
CMB. For the first time, we introduce parametrization for non-linearities of the instrument for both temperature and polarization 
systematics. We have shown that although linear systematics fields can distort existing primordial bispectrum, they do not 
generate new spurious bispectrum. Perhaps more important, we show that certain types of second order non-linear systematics 



15 



generate spurious bispectrum even if the primordial CMB is perfectly Gaussian. Quantitatively, we propagate the effect of 
systematics errors to the CMB bispectrum and then to eventually assess their impacts on the non-Gaussianity parameter / NL . 

We consider different bispectrum configurations which are predicted by various well-motivated inflation models: local, equi- 
lateral, and orthogonal shape of non-Gaussianity. For each type of non-Gaussianity, we calculate the effect of instrumental 
systematics on the temperature bispectrum TTT, and bispectrum involving polarization fields: EEE, TEE, and TTE. We dis- 
cuss that by optimally designing the scan strategy, the level of several polarization systematics can be reduced. However, the 
quadrupole leakage can not be averaged out even with perfect coverage of rotation angles. Therefore, for polarization systematics 
we primarily focus on non-linear quadrupole leakage q. 

We calculate the tolerance limits on the systematics parameters so that spurious f^T'' does not degrade our ability to con- 
strain primordial non-Gaussianity. We discussed the implication of our findings to future CMB experiments in particular the 
space-based PLANCK mission. We find that the non-linear systematics does not significantly affects CMB power spectrum 
measurement. The effective local f^' could reach (9(10) before systematic effects show up in the CMB dipole measurement. 
As a result, if the non-linear systematics is not controlled by dedicated calibration, the measured bispectrum could be biased. The 
effect of linear systematics on the primordial bispectrum is small. For linear systematics control of rms ~ 10%, measurements 
at ( < 2000 do not suffer significant degradation, however for experiments probing £ max > 2000, even linear systematics could 
alter the cosmological bispectrum. Since secondary anisotropics start to dominate at around I ~ 2000, we conclude that linear 
systematics is much smaller worry than the nonlinear systematics. 
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Appendix A: Discussion on Polarization Systematics 



In the paper we mostly focused on the Jones parametrization formalism and generalized the formalism presented in Ref [3j 
However for the PSB it is convenient to parametrize systematics starting from Stokes parameters J48ll /, Q and U ( V — for CMB 
polarization) as opposed to the Jones matrix formalism which is particularly useful for describing coherent polarimeters ll32ll49ll . 
Systematics parameters in this formalism are: the gain factor g, the differential beam-width of the beams fi, the differential 
pointing p, the beam ellipticity e and the beam rotation e. These parameters are derived from differences in intensity in the 
Gaussian 2-D polarized beam response function for each polarization^ Although throughout this paper we employed the Jones 
formalism, it is constructive to compare it to the Stokes formalism 13911 formalism. This will shed light on the similarities and 
differences between the two. Table III summarizes the correspondence between the parameters of two formalism. 

Our starting point is the Fourier-space analysis of JUt. We expand the leading terms to first order and identify each one of 
them with the corresponding term of the Jones Formalism ll32ll . Since in real space the temperature and polarization patterns 
are convolved with the beams, these expressions are simply the product of their Fourier transforms in Fourier space. We restrict 
the discussion to an elliptical Gaussian beam (with major and minor axes <x A and cr y ) 



B(x) = 



1 



2no~ x o~ Y 



exp 



( (x- Px ) 2 (y-py) 



2\ 



and its Fourier transform is 



2o-2 



(Al) 



B(l) = exp 



+ A ■ p 



(A2) 



The pointing error merely shifts the phase of the beam representation in Fourier-space. 
From Eq. (12) of l39ll the beam function with ellipticity e and pointing p in /-space reads 



CO CO 



B(Y) = e-y J] Yj i2m+ni >n{z)J n {lp) x e 



i(2m+n)tjj-in8 gi{2m+ri){<pi—a) 



(A3) 
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Angles Layout 




N 

FIG. 6: Angles layout: The X-Y system is fixed to the focal plane, which is itself rotating with an angle /? with respect to the fixed N-S system 
on the sky. The X'-Y' frame coincides with the ellipse principal axes and the angle if/ represents the angle that the polarization-sensitive-axis 
makes with the ellipse major axes. In the presence of a half wave plate (HWP) the angle if/ varies, otherwise it should be fixed. The pointing is 
represented by the 2D vector p which makes an angle 6 with the positive X-axis. The angles represents the tilt of the ellipse major axis with 
respect to p. 



where 

f 

y = 4^ + ^) 

and e = ^' + ^ v and the pointing is defined by p x = p cos 6 and p v = p sin 8. The angle a is scanning strategy- and time-dependent 
and is shown in Figure 6. 

The observed CMB polarization obtained as in Eq. (18) of j39|] 

Q'±iU' = ^(B + )(Q±iU) + ^(B-e ±2ia )T+^(B + e ±4ia )(Q + iU) 

- ^{B + )(e ±4k ')(Q + iU) - ^{B_e ±2ia ){e ±4ia )T - L(B + e* 4ia )(e ±4ia )(Q ± iU) (A5) 

where angular brackets (...) represent averaging over 'hits' at the given sky-pixel. A scanning strategy is said 'ideal' if (e ±ma ) = 
for any n + and is 'uniform' if {e ±ma ) does not vanish, yet it is constant across the sky. The symbols B ± stand for \{B\ ± B2) 
where B\ and B2 are the first and second beams in the beam pair. The map inversion becomes singular when D = 1 - {e 4,a ){e~ 4,a } 
vanishes. This implies that a good scanning strategy should minimize (e 4la ) and (e~ 4,a ), terms responsible spin-flip, i.e. Q+iU — > 
Q ± iU, since (e 4la ) are spin-±4 fields. We therefore neglect the second line of Eq.( IA5b , noting that these terms need not vanish 
in general and approximate D » 1 to first order in scanning strategy quantities. Note also that we ignore pixel-rotation for the 
sake of simplicity since at the leading order it decouples from other beam systematics and can be accounted for by simply setting 

Q! ± iU' = e ±2ie (Q ± iU). (A6) 

The analog of Eq.(|A5]l in the |32| formalism reads (Eqs. 17 & 18 of JH) 



Q' + iU' = Q + iU + cr(p ■ V)(g + iU) + <r{d\ + id 2 )(di ± id 2 )T 

+ cr 2 q{dy ± id 2 ) 2 T + (a + 2iaS){Q ± iU) + {fx + if 2 )(Q + iU) + ( 7l + iy 2 )T. 



(A7) 
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Gain 
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Rotation 


CO 


Pixel rotation 




Pointing 


Pa 


Pointing 


£<cos(«)) + g:<sin(a)> 
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Pointing 


Pb 


Pointing 


£<sin(a)> - £<cos(a)> 


Flip 


fa 


Flip 


(cos 4o-> 


Flip 


fb 


Flip 


(sin4a> 


Monopole 


Ya 


Diff. gain 


g(cos 2a) 


Monopole 


Yb 


Diff. gain 


g(sin 20-) 


Dipole 


d a 


Diff. pointing 


^<<?'"> 


Dipole 


d b 


Diff. pointing 




Quadrupole 


q 


Diff. ellipticity 


e ■ e 2i * 



TABLE III: Comparison between the l32ll and 13^1 formalisms: the various parameters used in the two formalisms are identified. It is shown 
that most of beam systematics vanish in case of ideal scanning strategy (pointing, spin-flip, monopole, dipole). The respective names and 
parameters are those used in the l32ll and l39ll papers. 

Note, in particular, that since T is a spin-0 field it must couple to spin±2 to generate Q ± ill'. Now, since d\ ± idi is a spin±l, 
d\ + idi must be a spin±l field which can come from a non-vanishing dipole-moment of the scanning strategy (e ±la ). Similarly, 
d\ + ido is a spin+1 field and therefore must couple to a spin±3 field, such as the octupole moment of the scanning strategy, to 
form the required spin±2 field for temperature leakage to polarization. Expanding Eq. flA5t at leading order we can identify the 
ll32ll with the corresponding l39ll terms. This is summarized in Table III. 

Note also that the scanning strategy functions, e.g. (cos 4a), can be space-dependent and very much like to CMB lensing 
this may induce higher-order correlations. However, since the scanning strategy is completely uncorrelated with the underlying 
temperature and polarization this will show up only at the trispectra level and is therefore irrelevant for the current work. As 
mentioned above, the main difference between the formalisms stems from the different treatment of the polarization direction 
iff + X w i m respect to the pointing (Fig. 6). A good example where this difference can show up is the inclusion of a HWP in 
the experimental setup. Averaging over the angle 0_by employing a fast rotating HWP will result in no polarization in the ll39ll 
formalism while the parametrization adopted by Il32ll fixes iff +x to 0, i.e. the polarization sensitive directions are aligned with 
one of the ellipse principal axes. As has been shown by J39ll the level of systematic B-mode (for a given beam ellipticity) 
depends on iff + x and it maximizes for iff + x = 45°. In general, these two formalisms agree upon making the appropriate 
identifications of the model parameters. 

Although we have considered a pretty general class of instrumental systematics, there are several effects one may want to 
include to extend our analysis: leakage of galactic temperature signal to the polarization maps due to bandpass mismatch, beam 
side-lobes and variations in instrumental response with frequency. In our instrumental systematics model, we have assumed 
monochromatic response of the detector. Another concern may be that our toy systematic model does not capture all the realistic 
effects. For example, pointing errors would generally introduce polarization leakage regardless of the shape of the beam. A good 
model for the pointing error depends on the scanning strategy, as the pointing offset could be, in principle, strongly correlated 
for which the error keeps constant over a period of time. The pointing error would cause the misplacement of time-ordered-data 
samples into wrong pixels during map-making. In other words, the error increases the effective size of the pixels with extra 
smoothing effect (related to the rms of the pointing error), making every pixel on the produced CMB map have its own effective 
beam. 

If the bias of bispectrum introduced by instrumental systematics is noticeable, it is possible to correct for the instrumental 
systematic effects by estimating the bias induced by the systematics if we have a precise enough knowledge of them. One way to 
do this is to use the recovered temperature and E-mode maps as input sky and simulate the instrument using the known beams. 
The output of this simulation contains some non-Gaussianity, e.g. certain type of bispectrum, not present initially, which is an 
estimate of the observed spurious bispectrum. The efficiency of this correction obviously depends on how well the systematics 
are known and modeled. 



Appendix B: Full-sky Treatment 

In this section, we focus on the derivation of the systematics distorted/generated CMB bispectrum under the more appropriate 
spherical sky treatment. The full-sky formalism can be obtained by simply replacing the Fourier components with spherical 
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harmonic multipole moments. We follow the mathematical notation and structure presented in 150115 111 

* ® /m + j dhY',"*a(h)®(h) + J dhY', m b(h)@(n)®(h) 

jmm'm" . \ 1 j jmm'm" m'" r\ 

a I' m' 1 //' l" ~r / . °V"m'" J //' /" /"' \Jl" f m 



= 0, 



h Z Z 



P"m" 



note that the integrals over the spherical harmonics are replaced by the geometrical factors defined in the following 

jmm'm" _ I ji ym* ym' ym" 



jmm m in 
J ll'l"l"' 



= j dhY'f'YfYpYfi" . 
The CMB temperature bispectrum with systematics contaminations can then be expressed as 



nobs _ V" 1 

°hhh ~ Z_i 



l\ h h 
mi mi m 3 



/f\obs f\obs f\obs \ 



(Bl) 

(B2) 
(B3) 

(B4) 



where the average is over the CMB and systematics field realizations, leading to 



tfhh = Z \ h h h I [<®^,®^*W + Z Z Z <®^,®fcm J ®^3®/>^;r<>- / TS? mr + 2 Pem - 

y m m m-ljV v,mn,mm 33 3 



+ 2j2j2j2j <0/ ' m ' 0/ >2 0/ :>" fl « a «>4« 7 / 3 « ; + 2 Perm. 



(B5) 



Here the first term is the desired primordial bispectrum, second terms is the bispectrum due to non-linear response of the 
instrument and the last term is the distorted bispectrum due to linear systematics a. We should the details for the linear- 
systematics, but can similarly be applied to the non-linear response. Noting that the Wigner-3 j symbol obeys the identity 



y ( h h h\ 



l\ I2 I" 
mi ni2 m'{ , 



(2h + 1) 



we can re-write the distorted bispectrum as 



Kkk = B w 3 + Z B lvm Z C '" R + 2 Perm - 

i"i" 1 



where 



R 



- S I 



h h h 
mi m 2 my 



In order to evaluate R, we first re-express I"],"?," 1 " as 



HI' l" 



fin 



h 1' 2 ' l'l 
K mi m'2 



( i r i" 



jinnm^m-, r m}m 2 m 3 
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(B8) 



where 



Then the expression for R can be re-written as 



///'/" 
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where, in the last step, we have introduced the Wigner-6 j symbol. The values of the Wigner-6 j symbol can be computed 
numerically with a fast and efficient recursive algorithm. 

Finally, substituting the expressions for R and including all permutations in a single expression, we can write the systematics 
distorted bispectrum as 



B hhh ~ B khh + Zj C f 

Ipq 



fhlpfhlq(-l)" ■ 



h h h 

I q p 



B t Pt ,fhipfhiq(-iy- 



l\ I2 h 

p 1 1 



B& pl iq fhlpfhlq(- 1 )" 



h h h 
q p I 



B 



pqh 
(B12) 



where n = (I + p + q). 



